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We study the classical dynamics of a quasiperiodic kicked rotor, whose quantum counterpart is known to be an 
equivalent of the 3D Anderson model. Using this correspondence allowed for a recent experimental observation 
of the Anderson transition with atomic matter waves. In such a context, it is particularly important to assert the 
chaotic character of the classical dynamics of this system. We show here that it is a 3D anisotropic diffusion. 
Our simple analytical predictions for the associated diffusion tensor are found in good agreement with the 
results of numerical simulations. 

1. Introduction 

One of the most remarkable interference effects is certainly the Anderson localisation Q): A 
classical random walker displays diffusive dynamics. But a quantum-coherent one (i.e. which keeps 
a memory of its phase) can either diffuse, or localise. This depends on the relative importance of 
diffusive transport and interference effects. 

Most studied wave systems which may show localisation are those where a wave propagates in a 
disordered medium, where disorder induces a random walk. However, the phenomenon of random 
walk is not restricted to disordered systems, as it can occur in very simple chaotic systems such 
as the well-known kicked rotor described by the Hamiltonian: 



where p is the momentum conjugated to the variable x (which, as it appears only in the argument 
of a cosinus, can be restricted to [0, 2tt]), the particle mass is taken equal to unity and K represents 
the amplitude of the sinusoidal potential. This system has one degree of freedom, but the presence 
of a temporal forcing and its non-linear character allow the existence of a chaotic regime. More 
specifically, for values of the so-called stochasticity parameter K > 4, the classical dynamics of 
that system along p is a random walk: the walker " particle p" diffuses 

In 1995, Raizen and coworkers realised a kicked rotor by exposing laser-cooled atoms to a 
pulsed standing wave Q), which is "seen" by the atoms as a sinusoidal potential acting on the 
centre of mass variables. However, for cold enough atoms, the de Broglie wavelength becomes 
comparable to the spatial period of the standing wave, in which case their dynamics is quantum 
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(wavelike) and one can observe the phenomenon of "dynamical localisation" (0, Q)> that is, the 
freezing of the diffusive transport of the walker in the momentum space. This phenomenon has 
been shown to be equivalent to the Anderson localisation in ID ([j-O)- An impressive surge of 
theoretical and experimental work on the subject followed (f^4^3). 

The correspondence between the kicked-rotor and the Anderson model can be extended to 
higher dimensions: It turns out that a quasiperiodic kicked rotor with d incommensurable tem- 



poral frequencies is equivalent to a d dimensional Anderson model (IffiHffa). Adding two new 
(incommensurate) frequencies to the Hamiltonian ([T]) thus prompts to the observation of the 
Anderson transition, which exists only in 3 (or higher) dimensions. We have engineered in our 



cold atom experiment (\24 - \2a ) a quasiperiodic realisation of the kicked rotor described by the 
Hamiltonian 

p 2 

H qp = — + K cos x [1 + ecos(u)2t) cos(uj3t)] S(t — n) . (2) 

n 

to study and characterise this transition, including, in particular, the first non-ambiguous meas- 
urement of its critical exponent. Here, e is the amplitude of the quasi-periodic modulation of the 
strength of the kicks, and the frequencies cj 2 and UJ3 must satisfy a condition of incommensurab- 
ility. 



It can be shown \2a . \23 . \2d ) that the above ID system ([2]) has the same transport properties 



as that of the following 3D kicked "rotor" 



P 2 \ ^ 

#3 = y + W2P2 + u 3 p 3 + Kcosxi[l + e cos x 2 cos x 3 ] S(t - n) , (3) 



where p = (pi,P2,P3) is the momentum conjugated to x = (x±, X2, x$). In fact, the temporal 
evolution of an arbitrary initial condition ip qp (6,t = 0) for H qp ([2]) is identical to that of the 
corresponding tp 3 : 

M0,t = 0) = Vqp(0i,i = 0)8(6 2 )5(e 3 ) , (4) 

for H 3 ([3|) (see |H)). Thus, our experience with the quasiperiodic kicked rotor can be regarded 
as an experience of transport of waves in a 3D medium originating from a time-pulse 6(t = 0) of 
a plane source (invariant along the directions P2 and ^3), where transport properties are actually 
observed in the direction perpendicular to the emitter. This presents strong similarities with the 



methods of study of transport of classical waves in 3D disordered media §21 ). 
Note that ([3]) is not the standard 3D kicked rotor (f7j ): 

222 

r Y) T) T) 

H 3 = -y- + ui2 + UJ3 + K cos x\ cos X2 cos X3 <5(t ~ n ) ■ (5) 

n 

Two characteristics of ([3]) are noteworthy: First, the linear dependence of the transverse kinetic 
energy in p 2 and p 3 . Second, the "true" spatial direction "1" and the transverse "virtual" directions 
"2" and "3" present an anisotropy controlled by the parameter e. How these features affect the 
classical dynamics of ([3]) as compared to that isotropically diffusive of ([5]) is the subject of this 
study. 

Our motivations are the following. The linear dispersion of the transverse degrees of freedom 
in the Hamiltonian ([3]) seems, at first glance, to imply that only the direction "1" in ([3]) is 
disordered while the transverse directions "2" and "3" are quasiperiodic (see (0, [IIS (13)). This 



quasiperiodic character could prevent the classical transport in the transverse directions, thus 
making the localisation observed in the quantum regime not a consequence of subtle interference 
effects but attributable to the peculiar classical dynamics of this system (25). This is however 
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not the case. The classical dynamics of ([3]) is, as we shall see, fully diffusive in all directions, and 
the localisation observed in the quantum regime is therefore a quantum effect of the same nature 
as Anderson localisation. Second, the anisotropy of the classical transport is known to greatly 
affect the properties of the localisation, including a non-negligible dependence of the transition 
threshold on the degree of anisotropy ( e.g on the parameter e) (fg^ ) . We indeed observed such a 
dependence (see the phase diagram in (jga ))- 

2. Analysis of the diffusive dynamics 

The Hamiltonian ([3]) being periodic in time, its classical dynamics can be studied by adopting a 
stroboscopic point of view, e.g. by focusing on the system state just after the kicks. Integrating 
Hamilton's equations from a kick to the following one, we obtain the equivalent of the Standard 
Map U) for this system, which relates values of the space and momentum vectors (x n p n ) just 
after the kick to those (x n+ i,p n+ i) after a period: 

Pl„ +1 = p\ n + K sin#i n (l + ecosx2„ cosx3 n ) , 
P2 n+1 = P2 n + Ke cos x\ n sin X2 n cos x% 
P3„+i = P3 n + Ke cos x i n cos x 2 n sin X3 
xi n+1 = x hi +pi n+1 , 

X2 n+1 = X 2n + UJ 2 , 
X3 n+1 = X ln + UJ 3 . 

Starting from an initial condition (xq,Pq), one can easily simulate the evolution of this system. 
In particular, we want to see how a state initially localised at po = expands as time passes. 
Afterwards, only properties averaged over many initial conditions po = and Xq £ [0, 2n) will 
be considered. 

Note that the temporal dependency ([6]) of the transverse variables X2 and X3 is trivial, and 
very different from that obtained from ([5]): 

X2 n + 1 = X 2n + U2P2 , 
£3„ +1 = X ln + U 3 p 3 . 

However, from © we can express J>2 n +i as: 

n 

P2„+i = P20 + Ke ^2 cos x h sin(x 2o + cos(x 3o + uj 3 j) , (8) 
i=i 

which shows that the chaotic character of the motion in direction "1" progressively "diffuses" into 
the transverse directions, which hints for a complex dynamics in all directions. Indeed, as in the 
case of the ID kicked rotor, the variables X\ j can be considered as uncorrelated random variables, 
uniformly distributed over [0, 2ir) (the quasi- linear approximation is valid for large K) and p 2 
thus also performs a random walk 

The randomisation of x\ is due to two conjugated phenomena: the sensitivity on initial condi- 
tions and the folding of the phase-space. The sensitivity on initial conditions is characterised by 
the Kolmogorov-Sina'i entropy h (an average of individual Lyapounov exponents over a chaotic 
part or all of the phase space - see (0)), which gives the time-scale 1/h over which two trajectories 
initially close have diverged by 5x\ ~ 2ir. The quasi-periodic action of X2 and X3 is to increase 
the entropy h, as compared to the standard case of the kicked rotor. Thus, we can say that the 



(6) 



(7) 
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variable x\ is at least as random and uncorrelated as the variable x of the kicked rotor (jgil ). 
Indeed, we know that, for K > 4, the whole phase-space of the kicked rotor is chaotic and that 
x can be considered as an uncorrelated (or shortly correlated, the correlation time being ~ 1/h) 
variable. This will be true also for ([6]). However, for e not too small, the threshold in K should 
be smaller. We have not dealt with this issue of the threshold of the diffusive transport in the 
3D kicked "rotor" ([3|). In the following, we will restrict ourselves to values of K larger than 4. 

The dynamics of the system in momentum space is thus diffusive, characterised by a diffusion 
tensor Dij = 1,2,3) defined as: 

Dij = lim^ (9) 
n 

where (X) represents the average of X over initial conditions. We can evaluate approximately this 
diffusion tensor. Let us consider the behaviour of (p2 2 ), f° r example. We can write this quantity 
as: 

n-l 

(P2 n 2 ) = K 2 e 2 ^2 {cosxijcosx lm ) xio x (sin(x 2o + u 2 j) sm(x 2o + u 2 m)) X20 x 
j,m=l 

x (cos(x 3o + u 3 j) cos(x 3o + U3m)} X30 . (10) 

The correlation function (cos x±j cos x\ m ) Xl0 decreases exponentially in \j — m\, as in the case 
of the ID kicked rotor So the main contribution to the sum in (|10p is given by j = m. 

Performing the average over the initial conditions, we finally obtain: 

K 2 e 2 

(P2n 2 ) ~ -g- X n , (11) 

in the quasi- linear approximation. Following those lines in the case of the other quantities (piPj), 
we see easily that the diffusion tensor is diagonal and reads: 

(12) 
(13) 
(14) 
(15) 

From this, it is clear that the diffusive transport is anisotropic and that the axes "1", "2" and "3" 
are the principal axes. 





a (K 2 /2)(l+e 2 /4) 


D 2 2 * 


a K 2 e 2 /8 , 


D33? 


a K 2 e 2 /8 , 




s . 



3. Confrontation with the results of numerical simulations 

In this section we will compare the predictions of the previous sections with numerical simulations. 
First, we verify that the transport is indeed diffusive in all directions (for K > 4), and that the 
non-diagonal coefficients D^j of the diffusion tensor vanish. Figure [1] represents the temporal 
evolution of (piPj) for the Hamiltonian ([3]) with UJ2 = lityfh and 0J3 = 27T\/l3 incommensurate 
with each other and with 2ir (this condition of incommensurability corresponds to the case of the 
quasiperiodic kicked rotor), K = 10 and e = 0.8. The time evolution of the diagonal terms (j>i 2 ) 
is linear while the non-diagonal terms (piPj^i) stay very close to zero. Also, the distribution of 
each pi has a Gaussian shape, characteristic of a diffusive transport (see figure [2]). 
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Figure 1. Time evolution of the correlation of classical momenta (piPj). The non-diagonal terms i ^ j (orange, blue and 
maroon) keep small values all along the evolution, whereas the diagonal terms increase linearly with time, the slope of (pi 2 ) 
(represented in black) being much larger than the slopes of (p2 2 ) (in red) and (p3 2 ) (in green). This shows that the dynamics 
is anisotropically diffusive and that "1", "2" and "3" are its principal axes. The dynamics is approximately the same along 
the axes "2" and "3" (see figure 0. The parameters are: K = 10, e = 0.8, ui2 = 2tt^/E and uig, = 27r\/l3. 
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Figure 2. Final momentum distributions along pi (in black) and P2 (in red), the distribution along p$ being approximately 
identical to that along p2- After 1000 kicks, they display all a Gaussian shape characteristic of a diffusive motion. The blue 
and green curves are fits by a Gaussian which do not show any statistically significant deviation. Parameters are K = 10, 
e = 0.8, ui2 = 2tt\/~E and 013 = 2n\/T3. 
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Figure 3. Dependence on K of the longitudinal diffusion coefficient Dn (black points) when the anisotropy is strong: 
e = 0.1. One observes oscillations around the average behaviour equation ^12) (represented in blue) which follow the known 
corrections of the kicked rotor, equation 11161 (dashed magenta curve), but only for values of K not too large. When K is 
increased, these oscillations decrease faster than in the case of the kicked rotor. The parameters are the following: u>2 = 2ttv5 
and LJ3 = 271-^/13. 




The dependence of the diffusion tensor versus K is now studied for e = 0.1, when the anisotropy 
is strong, and for e = 0.8, when the anisotropy is much weaker. 

For e = 0.1, Dn(K) is found to oscillate around its average behaviour, equation (112j) . the 
amplitude of these oscillations decreasing as K increases. This recalls the oscillations of the 
diffusion coefficient of the ID kicked rotor when K is varied. In the latter case, these oscillations 
appear due to the presence of residual correlations (sinx n sinxo). Taking into account correlations 
up to four periods, it leads to the following approximation for the oscillations of D(K) for the 
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Figure 4. Dependence on K of the transverse diffusion coefficients D22 (red points) and D33 (green points) when the 
anisotropy is strong: e = 0.1. The diffusion coefficients are found to oscillate around their average behaviour, equations H13H 
and Jl4j. The oscillating corrections do not decrease in amplitude as K increases, and they have the same period of the 
oscillations of D\\(K). The parameters are: L02 = 2tt\/5 and W3 = 27rVT3. 
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Figure 5. Comparison between the behaviour of D22 vs K for the Hamiltonian l(3j (red points - see figure [4} and for the 
Hamiltonian J5]l where X2 and X3 follow equations J7} (orange points). In the latter case, the oscillating corrections, which 
are a signature of residual correlations between kicks, decrease in amplitude as K increases. The parameters are: e = 0.1, 
ui2 = 2tt \/S and W3 = 2n\/T3. 



kicked rotor \3a\: 



D(K) « ^{1 - 2J 2 (K)[1 - J 2 (K)}} , (16) 

where J 2 {K) is the Bessel function of second order. In the case of the Hamiltonian ([3]), we expect 
a similar phenomenon to occur when e is small. D\i(K) should follow the same oscillations of 
D(K) (|16p . This is indeed the case, but only for not too large values of K, as can be seen in 
figure [3J As K increases further, the oscillations of D\\ (K) are quickly killed and we recover the 
average behaviour equation (I12p . 

The evolutions of the transverse diffusion coefficients D 22 and D33 vs K for e = 0.1 are 
represented in figure HI The oscillating corrections to the average behaviour equations (fT3l) and 
(|14p do not seem to decrease in amplitude. On the contrary, in the case of a true 3D anisotropic 
kicked rotor, these oscillations are found to decrease in amplitude (see figure [5]). The trivial 
dependence of x 2 and X3 in the case of the Hamiltonian ([3]) is certainly responsible for this. 

When the anisotropy is weak, e = 0.8, the residual correlations between kicks, responsible for 
the oscillating corrections on the diffusion tensor discussed above, are killed even for small K. 
This is clearly seen in figures [6] and [7] where the diffusion tensor elements are found to closely 
follow their average behaviours, equations (I12p . (I13p and (j!4j) . 

In conclusion of this section, we have shown that the classical dynamics of the model Hamilto- 
nian ([3]) truly corresponds to a 3D chaotic diffusion, anisotropic along the true spatial direction 
"1" and the "virtual" directions "2" and "3". The comparison of numerical simulations with the 
calculations of sect. [2] shows a very good agreement. 
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Figure 6. Dependence on K of D\\ (black points) when the anisotropy is small: e = 0.8. D\\ closely follows the average 
behaviour prediction, equation 11121 . represented in blue. In this case, correlations between kicks are all negligible. The 
parameters are: u>2 = 2n\/5 and ui-j, = 27rvl3. 



200 




Figure 7. Dependence on K of D22 (red points) and D33 (green points) when the anisotropy is small: e = 0.8. The diffusion 
coefficients closely follow the average behaviour predictions, equations H13H and 1141 . represented in blue. The parameters 
are: UJ2 = 2ir\/5 and ujg = 2-7r\/l3. 



Conclusion 



In this paper, we have reported some important statistical properties of the classical transport 
in the 3D system corresponding to the qu asip eriodic kicked rotor we have recently used to study 
the Anderson metal-insulator transition ([j^Hl)- The classical dynamics has been shown to be 
fully diffusive in all directions, a prerequisite of the possibility of 3D Anderson localisation in the 
quantum regime. We have given a simple analysis of the anisotropic character of the transport, 
which is found to be in good agreement with numerical simulations. This should enable us to 
gain more insight into the metal-insulator phase diagram of the quasiperiodic kicked rotor (fgjl ) . 
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